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Abstract 

A new fit is done to obtain numerical values for the order low-energy-constants 
L[ in Chiral Perturbation Theory. This includes both new data and new calculated 
observables. We take into account masses, decay constants, -ftr^4, tttt and ttK scat- 
tering lengths and slopes and the slope of the pion scalar formfactor. We compare 
in detail where the changes w.r.t. to the 10 year old "fit 10" come from. We discuss 
several scenarios for estimating the order constants and search for possible 
values of them that provide a good convergence for the ChPT series. We present two 
such sets. One big change is that the fits do not have the expected behaviour in the 
limit of large Nc as well as before. 



1 Introduction 



Since its very beginning Chiral Perturbation Theory (ChPT)[Tl El [3], the effective field 
theory of QCD at low energies, has been successful in the description of several hadronic 
observables. Unfortunately when one tries to perform loop calculations to improve the 
precision of the predictions, one faces a problem. The couplings appearing in the £4 
Lagrangian are many, i.e. 10, and they must be determined from phenomenology. One 
of the first determinations of such couplings was done already in [3]. There most of the 
next-to-leading order (NLO) couplings were inferred both from phenomenology and from 
considerations lead by large estimates (where N^. is the number of colours). 

Considering such good results it is important to decide whether ChPT is a suitable 
theory to achieve precise determinations of the hadronic observables. It urged then to 
carry on a program and perform next-to-next-to- leading order (NNLO) calculations [4]. In 
the last 10 years many two-loop calculations in three-flavour ChPT have been done, see 
[5] for a review. 

Notice however that going to higher orders raises a serious issue: the number of the 
unknown couplings increases rapidly. If on the one side adding loop diagrams should allow 
us to include better corrections and improve our descriptions, on the other side, many 
unknown parameters contribute and this seriously threatens the predictivity of the theory. 
Furthermore without knowing the values of such constants the convergence of the chiral 
expansions are difficult to test, although feasible with the method described in |6]. 

The two-loop expressions now available can be used to perform a new global fit at 
NNLO of some of the next-to-leading-order (NLO) low-energy-constants (LECs) into the 
game, the L[. A first attempt was done [7] when some experimental information was 
available for the Ki^ decay and by estimating the NNLO contributions using dispersive 
analysis. The fit was refined later on, when the full NNLO calculation for this process 
was performed jSJ I9|. After that many other observables have been calculated at NNLO 
and many of them are also better known experimentally. Therefore the time has come to 
perform a new fit of the L[ couplings at a NNLO precision. In this paper we present results 
for such a fit. Some studies using the extra observables but without performing a full new 
fit were reported in [TOl [HI [I2] • Notice that in a preliminary phase we converted all the 
FORTRAN programs used to evaluate the amplitudes up to NNLO into C++ code. Our 
fits are all performed using such programs. 

The paper is organized as follows. In Section [2] we sketch out the ChPT formalism and 
its main underlying ideas. In Section [3] we review the phenomenological input we used in 
our results. We also show which are the L[ that give the largest contributions for each 
observable we included. Notice that now much more input is present compared to the past 
fits [HI [9]. In Section m we present the main model we used to estimate the C[, i.e. the 
coupling constants appearing in the NNLO Lagrangian. Such an estimate is usually called 
the resonance estimate. In Section Owe summarize the status of the main existing NNLO 
fit so far: fit 10 of [9j. In Section E] we show our main findings using the estimates of 
Section HI We quote different fits so to show how each observable we include affects our 
findings. The best fit we get we call fit All and should be considered the main output of 
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this work. This fit exhibits several differences with fit 10. One especially striking feature is 
that it does not respect any longer the large Nc relation 2L\ ^ Lg. We also show for fit All 
the convergence of the expansions for masses and decay constants, which is much improved 
compared to the one of fit 10. In Section 16.11 we perform a fit of the L[ using as input 
different experimental results for the Ki^ amplitude. We show that with this input the 
large Nc relation 2L'[ ^ is better satisfied, although the resulting fit is not as good as 
fit All in convergence for the masses. In Section we try to justify and test our estimate 
of the C[. In Section [6^ we compare further the fits obtained using their predictions for 
the two-flavour LECs ii. Fit All again results as the most convincing one. In Section 1^3] 
we show results for fits based on a different estimate of the C[ couplings that can be found 
in [13], this is essentially a chiral quark model estimate (CQM). The fits are not as good 
as fit All, nor for the nor for the convergence of the expansions. Also in this case the 
large Nc relation is not satisfied. In Section 16.51 we provide a short comparison with the 
recent determination of of [14]. 

In Section[7]we show results for an another treatment of the C[. We let the C[ couplings 
to take random values, although they are forced to keep the size l/(167r^)^. These fits 
have been done requiring extra constraints of convergence for mass and decay constant 
expansions, as explained in Section 17.11 In this way it is easier to select only credible fits. 
The results of such a study are finally shown in Section 17. 2[ We found very many good 
fits that correctly predict all the observables used as input and with low x^- These fits 
are unfortunately different looking from each other. Therefore we can only show which are 
the ranges where we found the LI to vary. For some of the NLO constants such ranges 
are very wide. This method shows however that it is possible to fit the NNLO expressions 
to the observables with of the expected size and it also allows to study well the strong 
correlations between the couplings. Finally in the appendix we present a table where we 
quote our estimates for the NNLO couplings. 

2 Chiral Perturbation Theory 

We devote this section to a brief description of the formalism of three-flavour ChPT [Il[2l|3]. 
Introductory references are |15l [16]. The notation in the following is the same as in |1]. 
ChPT relies on the assumption that the flavour symmetry of QCD is spontaneously broken 
to the diagonal subgroup, SU{3)l x SU{3)ii — )■ SU{3)v- According to the Goldstone 
theorem, 8 pseudo Goldstone bosons then arise. These are identified with the low lying 
pseudoscalar mesons and are organized in a unitary 3x3 matrix 




(1) 



where (j) is a hermitian 3x3 matrix: 



TX' 



■+ 



(j) = n 





76^ / 



(2) 
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The Lagrangian describing the low-momentum strong interactions of the hght mesons must 
be invariant under SU{3)l x SU{3)fi local transformations. The most general lowest order 
Lagrangian is 

A = ^((^X) + (X+)), (3) 

with 

X± = u^XU^ ± ux'u , 
X = 2Bo{s + ip). 

The fields p, = v^ — a^ and = v^ + a^ are the standard external scalar, pseudoscalar, 
left- and right-handed vector fields introduced by Gasser and Leutwyler [21 [3] . The con- 
stants Fq and Bq are instead the leading-order (LO) LECs. The notation (X) stands for 
trace over up, down and strange quark flavour. 

Starting from this Lagrangian we can then build up an effective field theory by including 
loop diagrams and higher order Lagrangians, where operators of higher dimensions are 
included. Their coupling constants are then counter-terms, i.e. their infinities absorb 
the UV divergences coming from the loop diagrams. In this way one obtains a theory 
renormalized order by order. Unfortunately going to higher orders the number of operators 
allowed by the symmetries increase and therefore also the number of unknown coupling 
constants. At NLO there are 10 LECs, called L[, while at NNLO there are as many as 
94, called CJ". We will always quote the renormalized versions where the CI are made 
dimensionless by using the physical value of F^. The renormalization scale is chosen to 
he IX = 770 MeV. While there is in principle enough phenomenological information to fit 
the first ones, we still need to rely on theoretical models or on some other method for the 
latter ones, as those described in Sections HI 16.41 and [71 



3 Fitting procedure and input observables 

In this section we first show how we perform the fits and then we review shortly the 
observables we use as input and their values. 

3.1 

The fit is performed using MINUIT in its C-I--I- version [T71 [18] . MINUIT is a routine to 
find the minimum value of a multi-parameter function. The procedure to perform the fits 
is very similar to the one explained in The function to be minimized is the of the 
fit. It is obtained summing up the partial contributions Xi(part) QSich. input observables 

2 / •^'i(meas) •^i(cau\c) 

Xi(part) = I 
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j(part) (4) 

i 

where Xi(meas) are the physical values for each input observables and Axi their associated 
errors. Xj(caic) are the results as calculated by ChPT up to NNLO. The rest of this section 
is devoted to list the values and the uncertainties used for each Xj. 

3.2 Masses and decay constants 

The masses and decay constants of the light pseudoscalar mesons have been calculated 
at NNLO in [19]. We use them as physical parameters, namely as input to calculate the 
several observables, similarly to what was done in [S]. Their values are given in [20] and 
are 

= 139.57018 MeV, m^o = 134.9766 MeV, = 547.853 MeV, 
mK+ = 493.677 MeV, tuko = 497.614 MeV, (5) 
= 0.0922 ± 0.0002 GeV . (6) 

The measurements in ([5]) and (E]) differ slightly from the ones used in the latest full fit [9]. 

= 139.56995 MeV, m^o = 134.9764 MeV, = 547.30 MeV, 
mK+ = 493.677 MeV, m^o = 497.672 MeV, (7) 
= 0.0924 GeV. (8) 

Notice that preliminary results of our work have been reported in the proceedings [211 122] 
and several unpublished talks. Those results were based on the masses of [23], that differ 
slightly from both ©, ([6]) and ([7]), (jH]). 

All the new fits shown in this paper have been produced using the values in (l5|) and 
dHj). However we have not observed any substantial modification of the outputs when using 
the old masses and pion decay constant ([H]). As discussed below other changes in 
experimental input are behind the changes of central values. 

The masses depend at LO on Bom (with rh = {rriu + md)/2) and on Bonis. L^, Lg, Lg 
and Lg appear in the expression at NLO. In there is also a NLO contribution from L^. 
The decay constant depends instead on Fq, as an overall factor, and on and Lg. 

3.3 Fk/Ft^ 

As input observable for our fits we will use the ratio F^/F^^ to eliminate the dependence 
on the unknown constant Fq, since it contributes as an overall factor for F^ as well. In 
the end the value of then determines the value for Fq. 
The ratio takes the value [20] 

-i^ = 1.197 ±0.007, (9) 
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that is also in full agreement with several lattice estimates as reported in [23]. Using Fk/F^^ 
at NLO we are sensitive to Lg. To perform the fit we expand the ratio as 



Fk 
F^ 



1 + 



El. 

Fo 



+ 



El. 

Fo 



Fk 
Fo 



El 

ni Fq 



(10) 



NLO 



NNLO 



SO that we can keep track of the exact contributions from the different orders. However we 
also check that the same quantity estimated as 



Fq + Fk\p4 + FxlpG 

Fo + -^7r|p4 + -Fvrlpe 



(11) 



gives approximately the same valuqj. 

Notice the experimental result for the ratio Fk/Fj^ ([9]) differs substantially from the 
one used in [9]. Ref. [9] used Fk/ Ft, = 1.22 ± 0.01. The change in Fk/F^, is one of the 
major sources of difference with [Q] as will be shown later in Section [61 



3.4 The quark-mass ratio nis/rh 

For the masses we have a similar problem, they depend on the quark masses and on Bq. 
We thus use as was done in [H [9] the ratio of the strange quark mass over the isospin 
doublet quark mass m as input observable. The two following relations involving the light 
pseudoscalar meson masses hold at LO in ChPT 



rh 



m, 



Ott 



m, 



On 



m 



m, 



On 



(12) 



where with rriQ we indicate the meson masses at LO. They are calculated subtracting from 
the physical values the NLO and NNLO expressions. We include both relations in f ll2p in 
the fits. For the pion mass we use the neutral pion mass m^o. In the kaon case we need to 
correct the physical value for the mass since its electromagnetic contribution is sizeable. 
We take the average between mx+ and mj^o and then we subtract the electromagnetic 
contribution as stated by the Dashen's theorem and an estimate of its violation: 



m'xa.v 



[m 



K+ 



+ m^o - 1.8(m2+ - m^o)) = (494.50 MeV)^ 



(13) 



The factor 1.8 in ( 11311 is due to the corrections to Dashen's theorem where we use the value 
of [25]. 

The value of the quark mass ratio has been calculated by several lattice collaborations. 
The authors of [9] used as standard input nis/rh = 24 with a 10% uncertainty, but they also 



^We thank Veronique Bernard and Emilie Passemar for pointing out that these were significantly 
different for some of our prehminary fits. 
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checked that ms/rh = 26 was compatible. Here instead we use rris/rh = 27.8 as obtained 
by the Flavianet Lattice Averaging Group in [23], and we again adopt a 10% uncertainty 
for the error to be used in the fits when comparing with the theoretical values ( !T2l) . 

In the end the calculated NLO and NNLO masses are used to determine the lowest 
order mass or alternatively Botji. 

3.5 K£4 for mf actors 

The decay K^{p) 7r+(pi)7r^(p2)e^(p^)z^(pi/) is given by the amplitude j7] 

T = ^v:MPuh,ii - iMPiW^' - An. (14) 

In ( !T^ and can be parametrized in terms of four formfactors: F, G, H and R. 
However the i?-formfactor is negligible in decays with an electron in the final state. Using 
partial wave expansion and neglecting d wave terms one obtains for the F, G and the H 
formfactors [26j: 

F{s^, se, cose) = fsis^, si)e'^' + fpc'^" cose + ... , 
G{s^,Si,cos9) = gp{s^,se)e'^'' + . . . , 

H{sn,Si,cosd) = hp{s^,Si)e'^'' + . . . , (15) 



where we also assumed that the p phase is the same for the three formfactors. In ([T5|) 
s^(s^) is the invariant mass of dipion (dilepton) system, 6 is the angle of the pion in their 
rest frame w.r.t. the kaon momentum. The F and G formfactors were calculated at NNLO 
in [8]. The quantities are especially sensitive to L\, and Lg. Also the H formfactor 
is known at order [27j but we do not use it as input observable since it depends on a 
different set of LECs, those from the anomalous intrinsic parity sector. 

The measured observables are obtained by further parametrizing /^(sjr, si) and gpisj,, si) 

as 

gpis^^se) = gp + g'pQ^, (16) 

where = s-^/^imD — 1. f|T6l) can be used to fit the measured data points. In [9] the 
preliminary linear fit from the E865 measurement [28] was used as input. It has the values 

/, = 5.77 ±0.097, /; = 0.47 ±0.15, = 4.684 ± 0.092, = 0.54 ± 0.20. (17) 

Now more precise results from the NA48/2 experiment are available |29j and their second 
order fit of the formfactors read 

f f" a q' 

^ = 0.152±0.009, ^ = -0.073±0.009, ^ = 0.868±0.01, ^ = 0.089±0.02, (18) 

J s Is Js Js 
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Notice that in [29] no measure of fs is reported, therefore we always use the fs = 5.75±0.097 
from the E865 collaboration [30] in our fits. Multiplying by fs and combining the errors 
in quadrature we obtain as measures 

fs = 5.750 ± 0.097, fs = 0.874 ± 0.054, f'J = -0.420 ± 0.052, 

gp = 4.99 ±0.12, = 0.512 ± 0.121 . (19) 

In [29] there is also a linear fit which gives ^ = 0.073±0.004 and thus f, = 0.420±0.024. 
Deciding which of the two fits for the Fs formfactor should be used is a relevant issue. The 
problem is how much we can rely on the curvature of the formfactor F,. As a matter of 
fact it is difficult for NNLO ChPT to reproduce the large negative curvature as was 
also noted in [HI E]- A dispersive analysis approach combined to two loops ChPT, similar 
to the one done for tttt scattering [31], might clarify the situation. An indication of this is 
given by Figure 7 of [9]. It is visible there that the dispersive result for Ki^ decay [7J has 
a larger curvature than the two-loop result [9]. 

Let us conclude with a cautionary remark about the data. In [32j it was made clear 
that isospin breaking effects at threshold give important corrections. These have not been 
taken into account in the NA48/2 analysis [29], thus they might affect significantly their 
findings. 

3.6 TTTT scattering 

The TTTT scattering amplitude can be written as a function A{s,t,u) which is symmetric in 
t, u: 

A(7rV^ ^ TT V) = t, u) + 5''-''5^''^A{t, M, s) + t, s) , (20) 

where s, t, u are the usual Mandelstam variables. The three fiavour ChPT calculation of 
A{s,t,u) was done in [TT]. The isospin amplitudes T\s,t) (/ = 0, 1,2) are 

T°(s,t) = ?>A{s,t,u) + A{t,u,s) + A{u,s,t) , 
T^{s,t) = A{s,t,u) - A{u,s,t) , 

T\s,t) = A{t,u,s) + A{u,s,t) , (21) 
and are expanded in partial waves 

+ 00 

T\s,t) = 327i^{2i + l)Pe{cos e)tl{s), (22) 
e=o 

where t and u have been written as t = — |(s — 4m^)(l — cos6'), u = — |(s — 4m^)(l + cos6'). 
In f l22|) we indicate with Pe{cos6) the Legendre polynomials. Near threshold the are 
further expanded in terms of the threshold parameters 

tlis) = q'\a[ + h[q' + 0(g^)), = ^{.3 - Ami), (23) 
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"o 


0.220 ±0.010 


ml 

TV 


h° 


0.276 ±0.012 


_2 


«o 


-0.444 ± 0.020 


lO-^mO 


bl 


-0.803 ± 0.024 




a\ 


0.379 ±0.010 




b\ 


0.567 ±0.026 


IQ-^m-^ 



Table 1: The values of the scattering lengths and slopes as found in [31] and our fitting 
uncertainties. In the third column the normalization factors are given. We quote here only 
those scattering parameters added as input in our fits. 



where the scattering lengths, slopes,. . .. These thresholds parameters consti- 

tute our observables. Currently a very precise determination of these parameters exists. It 
is based on a dispersive analysis approach and on two-flavour ChPT and can be found in 
[31j . In Tabled] we quote the values of the threshold parameters we use in our fits and their 
corresponding uncertainties, which we took to be double the ones in [31]. For most fits we 
used only Oq and Oq but we have checked that the others listed in Table [1] are also well 
within the uncertainties quoted. Notice also that the NA48/2 experiment in [2S] obtained 
compatible values for and Oq from the measurement of the S = Sp — Ss phase shift in i^^4 
decays. 

3.7 ttK scattering 

The 7tK scattering process has amplitudes T^{s,t,u) in the isospin channels / = 1/2,3/2. 
They have been calculated at NNLO in ChPT in [12]. As for vrvr scattering, it is possible 
to define scattering lengths and slopes a^, h[. So we introduce the partial wave expansion 
of the isospin amplitudes 

T\s,t,u) = 167r^(2£± l)P^(cos^)t^(s), (24) 

1=0 

where P^(cos^) are the Legendre polynomials. Then we expand the tl[s) near threshold 

4{s) = l^Tsql^K (4 + bW.K + 0{ql^)) , (25) 

where 

= 4 1 ) V ^ ) ' ^ ^ 
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1/2 
«0 


0.224 ± 0.044 






0.85 ±0.08 


lO^^m;-'^ 


3/2 
«0 


-0.448 ±0.154 


lO-^m-i 


1,3/2 


-0.37 ±0.06 


lO-^m-^ 


1/2 


0.19 ±0.02 


10-^m;3 


hT 


0.18 ±0.04 


lO-^m^s 


3/2 


0.65 ±0.88 


lO^^m;^^ 


6f 


-0.92 ±0.34 


lO^^m"^ 



Table 2: The values of the scattering lengths and slopes as found in [33j. The uncertainties 
quoted here are those used in our fits and are double the ones of [33]. In the third column 
the normalization factors are given. We quote only those scattering parameters added as 
input in our fits. 



is the magnitude of the three-momentum in the center of mass system. The Mandelstam 
variables are given in terms of the scattering angle Q by 

t = -2g^^(l-cos^), M = -s-t±2m|^±2m2. (27) 

( |25|) defines the tiK scattering parameters a\ and h\ that are our input observables. These 
have been computed from Roy and Steiner type equations in [33]. The results for the s 
and p waves scattering parameters we use are reported in Table [2l For most numerical 
results we used only and aT but we have checked that all the others agree within 
uncertainties. 

3.8 Scalar for mf act or 

The scalar formfactor for the pion is defined as 

Fl{t) = {Ti\v)\uu^dd\T,\q)), (28) 
where t = p — q. Near t = it is expanded via 

= (0) (l ± lirYst + c-st' + . (29) 

The observables (r^)J and are used as input in our fits. The NNLO ChPT calculation 
for these quantities was performed in [10]. The scalar formfactor cannot be measured 
experimentally. Measuring the vrvr phase shifts and using a dispersive representation it is 
possible to infer its energy behaviour and therefore the values of (r^)^ and [Ml [5B] 

(r^)- = 0.61 ± 0.04 fm^, = 11 ± 2 GeV-^ (30) 
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Notice that the result for (r^)^ is also compatible with the lattice result of [37] . 

This is all the information we can extract from the scalar formfactors. Currently, there 
are basically no results available for Fg{0) and for the energy behaviour of the kaon scalar 
formfactors or of the strange contribution to the pion formfactor. 

3.9 LI and L\q 

We do not attempt to fit the remaining NLO LECs, Lg and L^q. Those LECs we fit are 
independent of Lg and L[q, or alternatively, none of the observable^ we discuss depend on 
them. One needs to include additional information to constrain their values. 

Lg appears alone at NLO in the electromagnetic radius of the pion vector formfactor. 
The NNLO contribution dependent on the other L[ is rather small [38j. It was therefore 
possible to fit that constant almost independently from the other couplings [38J . Further- 
more it never appears at NLO in any of the observables used here as input thus it does 
not affect much our fits. We always set Ug = (0.593 ± 0.43) x lO'^ for /x = 0.77 GeV. 

L^Q can be estimated using r decays data on the V — A spectral function [39j . Its value 
was found to be L^q — (~4.06 ± 0.39) x 10^^ at /i = 0.77 GeV. However this constant 
never appears in the observables under study not even at NNLO. Therefore it does not 
have any influence on our fits. For this reason we always set such constant to zero in our 
fits. 



4 Resonance estimates for the C[ 

The many unknown coupling constants that appear in the Lagrangian, the C^, represent 
the major problem for performing the fit with a 0{p^) precision. A lot of effort went into 
trying to estimate them using different models and treatments. The one we present here, 
also used in [S], is the resonance saturation model [lOl SI]. It is based on the idea that 
the LECs encode the information from physics above AchPT ~ 1 GeV, and that they are 
dominated by the physics just above this scale, i.e. the physics of low-lying resonances. 
Therefore we need a Lagrangian that describes these new particles and their interactions 
with the pseudoscalar mesons of the octet. We include only vector, scalar and the rj' fields. 
We use the same estimate described in , thus we refer the reader to that paper for further 
details, including the Lagrangians used at the resonance level. 

The model is used then to estimate the contributions depending on the C[. In [8], 
the heavier mesons were integrated out producing p^ Lagrangians for the pseudo-Goldstone 
boson. The heavy resonance fields for the vector mesons produce 

= -^^(VA(k,«1)K,x-]) + ^|^(K,/!:1(v^K,M)) 



^with the exception of K£4 where a very small dependence is present for S£ ^ 0. 
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and the scalar mesons 



((VA/r)K,x-]), (31) 



2 



(32) 



While for the t]' they obtained 



^v' = -^d,{x^)d'^{X-) (33) 
In ( 13T|) and ( 132|) /^'^ are defined as 

In [8] the above Lagrangians were not rewritten in the standard form of the Lagrangian at 
p^. That work has since been done using more general resonance lagrangians in [121 
We have checked that the results using the Lagrangians fl31f32p directly agrees with the 
same inputs using the CI' directly in terms of resonance parameters as derived in [l2l H3] . 
The T]' contribution was rewritten in the C]' in [H]. 

The values we choose for the different couplings are the same as in [8] 

fv = 0.20, = -0.025, gv = 0.09, 

av = -0.014, Cm = 42 MeV, q = 32 MeV, 

dm = 20 MeV. (34) 

and the masses are the experimental ones [2U] . 

mv = mp = 0.77 GeV, ms = 0.98 GeV, rUn' = 0.958 GeV, (35) 

In Table [TS] in the appendix we quote the C[ as estimated through the resonance model. 
We did not include more sophisticated resonance models because this would have again 
increased strongly the number of free parameters to be fitted. As discussed below we 
also have indications that terms suppressed by l/N^, the number of colours, might be 
important. These cannot at present be estimated using this type of approach. 



5 Existing fits 

In this section we describe a bit more in detail the earlier fits. The main full fit done is 
fit 10 in [9jfl. Earlier determination of the L[ did not fully include NNLO effects and we 
thus do not discuss them here. The values for the obtained in fit 10 are reproduced in 
Table [3] in the column labelled fit 10. This is a full NNLO fit of the L[ and it was done 
including the quantities and the L[ whose value they infiuence most: 

The E865 data were still preliminary then, the main fit in 9 was with older K14 data. 
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fit 10 


fit 10 iso 


10^ LI 


0.43 


0.39 ±0.12 


lO^Ll 


0.73 


0.73 ±0.12 


lO^Ll 


-2.35 


-2.34 ±0.37 


lO^Ll 


= 


= 


lO^Ll 


0.97 


0.97±0.11 


lO'Ll 


= 


= 




-0.31 


-0.30 ±0.15 


lO'Ll 


0.60 


0.60 ±0.20 


(dof) 


0.26 (1) 



Table 3: The results for fit 10 of [8] and for a similar fit done without including isospin 
breaking corrections for the masses (fitlO iso) and also using the masses in ([5]) and decay 
constant as in ([6]). The uncertainties are those calculated by MINUIT. The two fits 
reported are in agreement within uncertainties. 





p2 


4 

p 


p6 


ml 


0.753 


0.006 


0.241 


m\ 


0.702 


0.007 


0.291 


ml 
FjFo 


0.747 


-0.047 


0.300 


1 


0.136 


-0.075 


Fr/ Fq 


1 


0.307 


-0.003 


Fk/F^ 


1 


0.171 


0.049 



Table 4: The convergence of the expansion for the meson masses and the decay constants 
for fit 10 iso. A similar behaviour holds for fitlO. The masses quoted are normalized to 
the physical masses, while the decay constants to Fq (Fq = 0.0869 GeV). 
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1. masses and pion decay constant with the old values as in ([7]) and (|8]) 

2. the KfA formfactor parameters: fs, f'g, Qs, g's- They constrained mostly Ll,L2,L^. 

3. Fk/F^ = 1.22 ± 0.01, sensitive to L^. 

4. rris/rh = 24 constrains Lj, Lg via the masses in f|T2|) . 

5. L4 = Lg = since they are l/Nc suppressed couplings. 

The Ci contributions were estimated using resonance saturation as described in Sec- 
tion m They also included there the axial- vector resonances, although their contribution 
was rather small. The scale of saturation was set to /x = 0.77 GeV, but /i = 0.5, 1 GeV 
were within errors. In fit 10 isospin breaking corrections in the masses and decay constants 
were also included, though the authors of [9] noticed that the neglect of isospin violation 
was a good approximation. Indeed the fits performed including or not these effects are in 
agreement within errors as can be seen from Table E] comparing the columns fit 10 and fit 
10 iso. 

Fit 10 has been so far a quite successful fit. Not only because it already included many 
quantities at order p^, but also because the resulting L[ nicely confirmed the estimates from 
resonance models. These are lead by the large Nc expansion which predicts e.g. 2L^ ^ L2 
and LI ^ LI ^ 0. While the second relation was imposed, the first one was found to be 
well satisfied. This added credibility to the fit itself even though it relied on the resonance 
estimate for the tree-level contributions. 

However the convergence of the perturbative expansion for this fit is not as expected. 
The different orders for the masses and decay constants are reported in Table H] for fit 10 
isc0. The 0{p'^) order of the masses turns out to be tiny, far less than the expected 30%. 
On the other hand the NNLO contribution is definitely too large. The sources of this bad 
convergence are basically two. First the constraint L4 = Lg = that clearly sends to 
zero many contributions coming from the NLO tree-level diagrams. Secondly most of the 

appearing in the masses expressions are estimated to be zero as well. Therefore they 
cannot help in canceling large two-loop contributions. On the other hand the convergence 
for the decay constants is quite satisfying. 

After fit 10 was performed many other observables have been calculated at 0{p^) in 
SU (3) ChPT such as the tttt and nK scattering threshold parameters. Of course it is very 
important to compare the pure ChPT predictions obtained using fit 10 with the values 
of Tables [T] and [2J These comparisons have been done and can be found in Table 1 of 
^llj and in Table 4 of |12j. Fit 10 is mostly in agreement within errors, although there 
are small discrepancies in some of the threshold parameters. Some comparison with scalar 
formfactors was done in [10]. The last three papers used the same inputs as fit 10 and 
tried to vary L4 and Lg to see if some preferred regions could be found. Here we redo the 
fit from the beginning with all inputs. 

''The numbers for fit f itself can be found in Table 2 of [10]. 
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1 


1 


4 


4 
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Table 5: Several global fits compared to fitlO iso. For all the fit yU = 0.77 GeV. The errors 
quoted are the ones as calculated by MINUIT. The numbers in the row labelled dof are 
the degrees of freedom for the fit. See the description in the text for further details on how 
the fits have been performed. The column labelled All is our main new result. 



6 New Fits 

The aim of this section is to show how the new measurements and observables included in 
our global fits, change the results compared to fit 10. We have rewritten as mentioned above 
all programs into C±± and are using the isospin symmetric versions of the calculations. 
We therefore first redid the fit using the same inputs as fit 10. The outputs are given in 
Table |3] in the column labelled fit 10 iso. This also shows that the minor changes in masses 
and as well the isospin breaking corrections do not affect the fit values appreciably. We 
will now add the effects of the changed experimental inputs and of the additional inputs 
to see how they change the fitted values of the L[. 

In Table E] we present several fits. These have all been performed using the resonance 
estimate of the C[ of Section [Hand Table [18] in the appendix, setting the scale of saturation 
/i = 0.77 GeV. Furthermore, we used the new values of the masses and decay constant of 
(|5]) and ([6]). We remind the reader that the use of these new parameter- values affects the 
output only within the uncertainties. Hereafter we summarize the steps in which we have 
included the new information. 

NA48/2 The input observables and their values are the same as for fit 10 iso, but we use 
the new measurements in (fT9|) for the K^^^^ decay from the NA48/2 collaboration |29j . 
The new measurements lead immediately to a striking feature: the large N^. relation 
2L\ ^ L2 does not hold any longer. It even turns out that L2 < L\. Notice that, 
as explained in Section 13. 5[ the slope comes from a second order fit of the fs 
formfactor and therefore it differs from the one used in fit 10. In Section [6?T] we will 
present also results for the linear fit of the fs formfactor. 

Fk/ Same as fit NA48/2 but with the new value in ([9]) for Fk/Ft^. is mainly affected 
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and becomes smaller than in fit 10. As a consequence also the convergence of the 
decay constants expansion is worsened, e.g. F^/Fo|p4 ^ 0.134 while F^r/Folpe ^ 
-0.126. 

AIIt^t In this fit we include a few more observables, i.e. the vrvr scattering parameters af] 

1 3 

and 0,2, the nK scattering parameters Oq and Oq and the pion scalar radius (r)|. We 
also release the constraint = Lg = because now we have many more observables 
included. Unfortunately since none of the new observables involve information on 
the masses (and therefore on those two couplings), we still cannot achieve the precise 
values of L4 and Lg. They are highly correlated and MINUIT gives a very large 
uncertainty for those. Furthermore since now the contributions to the masses 

due to L4 and Lg are not zero, and Lg diminish. 

All This fit is very similar to fit AIIt^t. Here we adopt the new value for the quark mass ratio 
rris/rh = 27.8. However we find that using values for the quark mass ratio between 
27 and 29 does not change the results considerably. The constants L^, that appears 
in the rj mass, and Lg are strongly affected by this change. This is also relatively 
true for L4 and Lg. We also tried to perform the same fit but setting L4 = Lg = 0. 
The resulting fit is very similar to fit NA48/2 but it has a huge ix^ = 45). 

C[ = In the last column of the table we quote the fit obtained including the same input 
as for fit All, but setting all the = 0. This fit has been done to show how the 
different can affect the LI fit. Notice that the constants L\, and L3 change a 
lot, while the others stay in the same area as in fit All. This is not surprising: the last 
few constants are indeed primarly fitted from quantities where many contributing C[ 
are large- A^c suppressed and those which are not are set to zero also in the simple 
resonance estimate used. 

All Same fit as All but all expressions are now at NLO. Use this fit for one-loop ChPT 
results. Note that this produces very high values for L4 and Lg. The underlying 
reason is that the lower value of LV/L^ requires a smaller L5 than before and the 
pion scalar radius then requires at this order a larger L4. This effect is also visible 
in fit All but is reduced when including the NNLO corrections. 

Fit All is what we consider as the present best fit for NNLO ChPT calculations, it thus 
superseeds fit 10 of |9]. 

Let us discuss how the ChPT expansion is affected by the new values for the L^. The 
various terms of the mass expansions read for fit All 

0.084 m^lpG = +0.049, 
0.181 m^|p6 = +0.075, (36) 
0.224 mj|p6 = +0.038, 
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loops only 
p^-L1 fit 10 iso 
p^-Ll fit All 


-0.02276 -0.02287 
0.1141 0.01115 
0.00342 0.00330 
0.00232 0.00247 



Table 6: The results for /+(0) in the two Ki^^ decays. This is an update of Table 3 in 



while those for the decay constants are 
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(37) 



In ( !36l) and ( !37|) we used the same normalizations as in Table HI although now Fo = 0.065 
GeV, this is due to the larger value of L4 which comes however with a large error. Notice 
that the convergence of the mass expansions in ( l36l) is improved compared to the one of 
fit 10 in Table |H However (136|) looks strange: the LO masses are larger than the physical 
ones and there are significant cancellations between NLO and NNLO. Furthermore, even 
if the convergence is improved, it is still quite different from the one expected. E.g. the 
,4 contribution is much smaller than the expected 30% and it is of the same size as the 



m 



nip 



p order. The convergence for the decay constants is a bit worsened compared to the one 
of fit 10, due to the low value of L5, but it is still acceptable. Notice also that when the 
ratio Fk/F„ is calculated with ffTTj) the resulting value is 1.168, which is 3% smaller than 
the expected 1.197 . This can be due to higher order corrections that are included in the 
ratio of ([II]), but not in ffTU]) . 

We performed more fits than those quoted in Table O We included more tttt scattering 
parameters and tj'K scattering parameters. We found that these fits are compatible with 
fit All of Table [5] within uncertainties. The same is true when we add the quantity cj. 

For completeness we quote the value for /+(0), the formfactor of Kfs decay at zero 
momentum transfer. This quantity was calculated at NNLO in |45j. We check how much 
the new L[ of fit All would affect its value. Notice that at zero momentum the dependence 
on Lg drops out. The other L[ appear only at NNLO. We discuss both the case of the 
charged kaon, the K'^ decay, and the neutral one K^^. The results we obtain using the 
masses of the particles involved in the decay are given in Table [61 These numbers can be 
seen as an update of those in Table 3 of [45]. The numbers for the pure loops are changed 
w.r.t. |15] mainly because of the change in F^. Fit 10 iso is essentially the same as fit 10 
used in [15] but has a small (1%) difference. The L[-dependent contribution changes but 
since it stays small the total result is very similar to [15]. The total value for /+(0) for K^^ 
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is instead 

/+(0)lc,=o(Lr=fitio) = 0.9916, /+(0)|c,=o(Lr=fitAii) = 0-9908 (38) 
and for /+(0) for is instead 

/+(0)lc,=o(Lr=fitio) = 0.9921, /+(0)|c,^o(Lr=fitAU) = 0.9910. (39) 

By comparing the two results in ( 138|) and in ( 139|) one realizes that the differences due to 
the L[ have a very small effect on /+(0). 

The main uncertainty still remains the value of the contribution of the Cj. The estimate 
for the relevant constants used in this paper leads to /+(0)|^^ ~ —0.045 but the fitting 
inputs used here do not strongly constrain the relevant combination. 



6.1 Linear fit for Ki^ decays 

One of the most striking features of the results presented in Table |5] is that as soon as 
the new results from the quadratic fit of the NA48/2 collaboration p9] are included, the 
constants L[ and take unexpected values. Indeed, as was noted above, they do not 
respect any longer the large- iVc relation 2L\ ^ Lg, but already in fit NA48/2 they are 
while when also the vrvr and the irK scattering lengths are included (fit All) we 
even obtain > Lg. On the other hand when we calculate the curvature f'J using the 
L[ as obtained in fit All we obtain /" = —0.124 to be contrasted with the experimental 
value f" = —0.437. Furthermore, whenever we include as input also /" we again obtain fits 
compatible to the ones in Table |5l but with much larger (e-g- ~ 35 for fit All) and 
the largest contribution comes exactly from These results confirm what was already 
stated at the end of Section [3751 i.e. the state-of-art ChPT does not reproduce such a large 
negative bend. Since and /" are highly correlated, the linear and the quadratic fit of 
the Fs formfactor present rather different slopes. 

For such reasons we perform fits using the slope of the linear fit f'g/fs = 0.073 [29] 
as well. The resulting fit, analogous to fit All, is reported in Table [71 By inspection one 
can see that the large Nc relation 2L^ = L2 still does not precisely hold, but at least 
1.4L^ ^ Lg. On the other hand L4 and Lg are again not suppressed, while Ly and Lg are 
unexpectedly small. Also all the constants have large uncertainties. 

The convergence of the chiral expansions is worse than the one for fit All. The various 
terms of the mass expansions read 
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(40) 



while those for the decay constants are 
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-0.04 ±0.40 


W^Ll 


0.02 ±0.79 
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Table 7: The analogous of fit All, but the linear fit of the formfactors have been used. 
The values of the L[ are at = 0.77 GeV. 
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(41) 



In the light of these results it is rather difficult to draw a conclusion. The very different 
predictions for L\ and obtained in fit All and this fit confirm that the picture of ChPT 
for decays is still incomplete. As mentioned earlier, we expect a dispersive analysis to 
produce a larger curvature. 



6.2 Some small variations on fit All 

In the resonance estimate described in Section H] there is at least an assumption not entirely 
justified. We assume the scale at which the saturation happens to be 0.77 GeV, i.e. the 
mass of the lowest lying resonance. Nothing prevents us to choose a larger or smaller scale, 
although this is still expected to be in the same range of energy. To check whether this 
assumption is safe we try to fit from data the saturation scale parameter as well. The 
results are rather reassuring. Fit All of Table [5] is completely unaffected by this procedure. 
The fitted saturation scale is 0.77 ± 0.45 GeV. 

The fit in Table[7]shows a little difference. The fitted saturation scale is now 0.71 ±0.31 
GeV. However the L[ do not change that much and the look of the fit is pretty much the 
same as before. 

We also attempt to find better estimates of the CI constants releasing the values of the 
couplings gvi Cd and Cm- Again we try to fit them using both the input of fit All and of the 
fit in Table [7] (linear fit of NA48/2 instead of quadratic). In the first case we find in fact gv 
and c„ close to the ones in ([MD- They read gy = 0.097±0.123 and c„ = 0.045±0.049 GeV. 
For Cd we find instead a value larger than expected, i.e. q = 0.093 ± 0.100 GeV. Anyway 
they are all affected by large uncertainties. The L[ fit is somewhat compatible with the one 
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4.0(4.1) 


4.8(4.8) 
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4.4 ±0.4 



Table 8: The values of the scale-independent SU{2) LECs ii. In the first three columns we 
show the values as predicted by fit 10, fit All and fit All linear using the NNLO matching 
conditions of |l6]. The numbers between parenthesis are the NLO results. In the last 
column we quote the known values from [HT] . Notice that the uncertainty over £4 is 
double the one quoted in [31] due to the still unclear situation for the lattice results [24J. 

of Table O fit All, within uncertainities because of the large ranges allowed for gv,Cm,Cd. 
The fits are in a very broad minimum here with only one degree of freedom. 

If we apply the same procedure but with the same Ki^ input as for Table [7] (NA48/2 
linear fit) we arrive to similar conclusions: the values of gv and q are similar to the ones 
in ( 134|) . while Cd is larger. Again all the resonance couplings present large uncertainties. 
The L[ are here rather well compatible with those in Table [71 

We also try to multiply the C[ by an overall constant a and include it as a fitting 
parameter. It is encouraging to see that the result is a ~ 1.03, namely the best fit is 
reached with basically the same values of the C[ from resonance estimate. Obviously the 
fit obtained is very similar to fit All. When we apply the same procedure to the fit in 
Table [7] the constant a takes the value 0.90. This affects the fit of Table [3, but still within 
uncertanties. 

From this we conclude that fit All is stable against small changes in the resonance 
estimate of the C[. 

6.3 Adding input: £i constants in two-flavour ChPT 

The authors of [l6l HT] study three fiavour ChPT in the limit where the rris is assumed 
to be much larger than rh and the external momenta. In this case they can integrate out 
the strange quarks and ^^7(3) x SU{3) ChPT reduces to SU{2) x SU{2) ChPT. Matching 
the results from the two frameworks they calculate explicitely the dependence of the two- 
fiavour LECs (the scale-independent ii and the c[) on the strange quark mass and on the 
three-fiavour LECs. These relations have been worked out using two different methods at 
order in ^61 HT] . 

There exist different evaluations of the ii. £3 has been estimated rather well using 
lattice results [23]. ii, £2, ii and i^ have been obtained by matching two-flavour ChPT 
with dispersive results [31], but (contradictory) lattice results exist for those too [23]. We 
have increased the error on £4 because of this. 

In Table [S] we summarize all the values of these constants and the results obtained 
by plugging the L[ and C[ of flt 10 iso, flt All and flt All linear in both the NNLO and 
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NLO relations of ^6] W7\ . By comparison of the first three columns with the last one of 
Table [8] it is easy to see that none of the fits correctly reproduces all the ii values. A 
similar conclusion holds also for the fit in the next-to-last column of Table [SI where all the 
CJ' are zero. The disagreement with the ii in this case is actually even stronger. These fits 
encounter particular trouble in fitting £2- 

We tried to fit the ii whose values appear on the last column of Table [8] in addition 
to the inputs used for fit All. Not surprisingly it is not possible to accommodate all those 
inputs at the same time. The resulting is approximately 22 and its largest contributions 
come exactly from the ii. Excluding from the fit £2 but including the others improves the 
situation. The resulting L[ values are very close to the ones of fit All, the most important 
deviation being lO^Lg = —3.18. The takes the value 3.15 with 7 degrees of freedom. 
Also in this case the value for £2 is still far from the expected one. 

This is not surprising. In ^6] it was found that the constant £2 depends on the couplings 
L2, and on the combination 2Cl^ — Cl^. The authors of [46j observed there that to 
find agreement with the determined value of £2 the combination of C[ must not be zero. 
Unfortunately those are two large-Ai'c suppressed couplings and therefore they are set to 
zero in our resonance estimate (see Table [T8l) . We also try to fit those two C[ using also 
£2 as input observable, but this has been unsuccessful as well. In this way we manage 
to accommodate the value for £2, but then ii is off, since it contains also a different 
combination of C{^, and of Cg. This last coupling is also A^'c-suppressed and thus 
estimated to be zero. In the end there is no way out: when we try to fit Cg too, there are 
other quantities taking very different values. The C[ are too correlated to be able to fit 
only a few of them. 

As far as regards the fit in Table[7]the results for the ii are even less clear. Its predictions 
are reported in the third column of Table [HI It is straightforward to see that now even the 
predicted ii is off. Of course when we try to fit all the ii the is very large (x^ ~ 37.7). 
Contrary to what happened for fit All, the situation does not improve that much when we 
exclude £2- For this fit seems to be very hard to reach the correct value for ii too. The 
resulting in this second case is 5.82 with ii = —1.4. 

Finally an extra cautionary remark. Requiring that the SU (3) ChPT constants predict 
the values for the SU{2) ones might not be a very safe assumption. What it assumes is 
that both SU (2) and SU (3) ChPT work well for the same quantities. For the vrvr scattering 
quantities, which are very much determined by loop parts, relatively small differences can 
become amplified in the resulting values of the LECs. 

6.4 A chiral quark model estimate for the C[ 

We discussed above a simple resonance saturation estimate for the NNLO LECs C^. There 
are other attempts at predicting these values as well from chiral quark models. As a 
representative of this we choose [13]. It also is a large Nc approximation but with a 
somewhat different pattern than our resonance saturation. Their method is based on a 
study of the relation between the chiral Lagrangian up to order and QCD, they find 
as expected that the LECs can be given in terms of some Green functions of QCD. In 
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Table 9: The results as obtained using the CI estimates of [13]. Both the fits include the 
same observables as fit All of Table |5l In the second column the coefficient a multiplied 
by the C[ has been included as fitting parameter. The L[ are given at /i = 0.77 GeV. 
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Table 10: The results as obtained using the estimates of [13] showing the convergence 
for the fit where a is left free. The normalizations are the same as explained in Table HI 
Now Fo = 0.065 GeV. 

the evaluation of these Green functions, several assumptions and approximations are made 
such that it is not a full derivation but something like a chiral quark model. Their results 
are presented in Table IV of [13] . 

We also use their estimate to perform the fits. The results can be found in Table [9l 
There are results for two different fits. They have been obtained including all the observ- 
ables as for fit All of Table [51 In the first column we use the C[ as quoted in Table IV of 
[13] . whereas in the second column we multiply them by an overall constant a that is also 
fitted. This second fit was done because we observed that the values for the C[ of [I3] are 
somewhat larger than the ones of the resonance estimates of Table [181 The fit confirms 
this observation and finds as best value for a = 0.27, i.e. C[ considerably smaller than the 
ones in [13]. The value of the fo^^ the two fits of Table [9] is somewhat worse than for fit 
All. Indeed it seems that it is now very difficult to fit the slope g'p of the G formfactor for 
K£4^ decay. From the second fit of Table [9] one can notice that when the C[ are allowed to 



21 



6 



a 

X' (dof) 



Q [13] 
0.38 ±0.10 
0.88 ±0.12 

-3.20 ±0.47 
0.42 ±0.17 
1.62 ±0.77 
0.43 ±0.21 

-0.32 ±0.45 
0.92 ±1.07 

4.13 (4) 



Q X C[ [T3" 
0.35 ±0.11 
0.43 ±0.30 

-2.04 ±0.79 
0.91 ±0.35 
1.03 ±0.57 
0.87 ±0.77 

-0.11 ±0.34 
0.17 ±0.71 
0.22 ±0.47 
1.20 (3) 



Table 11: The results as obtained using the C[ estimates of [13]. Both the fits include the 
same observables as the fit in Table [TJ i.e. as fit All but with the linear fit from NA48/2. 
In the second column the coefficient a multiplied by the has been included as fitting 
parameter. The L[ are given at = 0.77 GeV. 
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Table 12: The convergence for the fit where a is left free is shown. The normalizations 
are the same as explained in Table HI Now Fq = 0.062 GeV. Fit as in rightmost column of 
Table M 

take smaller values the constant compensates for that. This allows the fit to reach a 
better value for the g'p. 

As can be seen in Table [THl even the convergence for the masses and decay constants 
is worse than the one for fit All reported in ( l36l) and ( 1371) respectively. Notice that we 
have not quoted the convergence for the fit obtained without multiplying the by the 
coefficient a. In fact this is found to be even worse than the one of Tabled the terms 
being constantly larger than, although comparable in size to, the ones. 

In Table [TT] and [12] we quote the results obtained with the C[ of [13], but fitting the 
slope of the linear fit for the Fs formfactor as in the fit of Table [71 Conclusions similar to 
the ones drawn for Table [9] hold here too. 



22 





m 


fit All 


best reso 


best rand 




0.76 ±0.09 
0.31 ±0.07 
1.10 ±0.14 


0.58 ±0.13 




1.40 ±0.009 
-0.99 
0.02 


1.40 ±0.009 
-1.06 
2.01 



Table 13: The values of the couplings Lg, Cl^ + C^g and Cl^ + '2.C{^ as obtained from a fit 
to lattice data p3] (second column) and from our fits. Notice that the resonance model 
used for fit All estimates the combinations of C[ occurring here as zero. The columns best 
reso and best rand are taken from Tables [TBI and [TSl and are discussed in Section [71 All the 
values are at the scale /i = 0.77 GeV. 



6.5 Comparison with a recent L\ determination 

The authors of propose a simplification of the NNLO predictions of ChPT to perform 
fits to lattice data points. As an example they fit the lattice results [Hj for Fi^jF^^ with 
their approximation to the two-loop ChPT prediction. Since too many couplings appear 
at NNLO they are forced to fit only a few of them and fix the others. They decide to fit 

and the only two C[-combinations contributing: C^^ ± and Cj^g ± 2C[y. C^g is 1/A^c 
suppressed and therefore set to in the resonance saturation model. C^^ and C\-j also do 
not get contributions from the resonances we included. In other models [H] C[4 and C\-j 
are estimated to be negative, but are very small in absolute value. In p3] the other L[, 
appearing at NNLO in Fk/Fj^, are set to the values of fit 10 [9]. The results of the fit 
to lattice data are quoted in the first column of Table [121 For comparison we quote in 
the same table the values as obtained by our best fits, including those for the random C[ 
search described in Section [71 and quoted in Table [181 in the appendix. The table shows 
that the value of in fit All is compatible with the result of [13] within uncertanties and 
even more compatible when we look at the fit Fk/Fj^ in Table [51 where we required = 
as in |T4:J. The fits best reso and best rand instead are very different. For these last two fits 
we can also compare the values of the C[-combinations. As far as regards the combination 

± Cis OUT random sets acquire both a negative value in contrast with f]3J. The second 
combination instead agrees upon the sign, but are rather different in value. 

For further comparison we can do also a fit similar to fit All but setting in addition 
lO^d^ = 0.31 and 10^C{^ = 0.55. The results are given in Table [H 

The fit obtained is different from fit All, in particular we find an even lower and 
the is somewhat larger. The value of Lg has decreased quite a bit while the others are 
compatible with fit All but the central values for L4 and Lg are rather large. Finally notice 
that the L[ of fit All are very different from the ones of fit 10 and might therefore affect 
the findings of [T^. We cannot draw any further conclusions at present. In future work 
we intend to include more lattice results which should clarify this issue. 
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m 


lO^Ll 
W^Ll 


n 4- n HQ 
0.53 ±0.21 

-2.97 ±0.43 
0.89 ±0.83 
0.30 ±0.09 
0.39 ±0.97 

-0.02 ±0.16 
0.13 ±0.20 


(dof) 


2.03 (4) 



Table 14: The fit as obtained using the same input as fit All and the as in the resonance 
estimate of Table [18] with the exception of and Cl^ set respectively to 0.31 x 10~^ and 
0.55 X 10~^ The L[ are given at /i = 0.77 GeV. 

7 Releasing the C[ 

All the fits presented in the previous section have unusual NNLO corrections to the masses 
and many also to the decay constants. In addition, if we included the requirement that 
the ii were also fitted well we could not find a simple good fit. 

An additional reason to go beyond what we have is that all the estimates used above 
with the exception of the singlet rj contribution only contribute to the NNLO LEC^ that 
are leading in N^. In the masses and decay constants in addition the estimates from 
the resonance exchange give no contribution at leading order in Nc at all. This is an 
unsatisfying situation, we do expect that the masses and decay constants should get some 
contribution from the NNLO constants. Inspection of the relations between the ii and the 
SU{3) LECs |16] shows that only combinations of the C[ appear that are suppressed by 
Nc- Thus especially the problem with £2 above requires some nonzero values for the 
suppressed constants. 

We could in principle allow all the C[ to be free and include them in the fit as well. 
However, from our earlier work in [B] it is clear that with the inputs used at present there 
are enough free combinations of the to fit all physical inputs directly. For this reason 
we also have explored another technique of estimate, based on a random wallj^ method. 
Hereafter we describe the main features of the algorithm used. See also the flowchart in 
Figure [TJ The algorithm is a version of simulated annealing. 

We first start with an initial set C[^'°'' = Cl^°^^\ These are chosen to be 

1. random (with a size given by l/3/(167r2)2 for those leading in l/N^ and 1/3 of that 
for the subleading ones), 

''This is true with the exception of terms involving (x-) which can get produced by the equations of 
motion. 

^The idea was born thanks to a discussion with Juerg Gasser and Gerhard Ecker. 
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^r(new) ^r(old) _^ ^r(stcp) 

and calculate x^^'"'^'' 




Figure 1: Algorithm used to select the random C[. It has been started with different values 
of the initial Cl^^'^\ as explained in the text. In the bottom decision square i? is a random 
number selected with a uniform distribution in the interval (0,1), while T is a parameter 
set such that it is of the same order of magnitude of the x^. More details can be found in 
the text. 
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2. all zero 



3. as obtained by resonance estimate (see Table fT8l) 

4. as obtained by multiplying the constants of [13j by 0.27 (see Table [18]). 

Then we perform the fit on the L[ using those c^.^°^'^\ After this we take a random step 
according to the formula 

^Knew) ^ ^Kold) ^ ^Kstep) ^ ^.(old) ^ _L_er, , (42) 

where in f H2|) is a random number generated through a uniform distribution in the 
interval (—1, 1) and we have used e = 0.01 and 0.001. For those that are Nc suppressed 
we further multiply d^'^^^^^ by 1/3. In this way the random generated set still respects 
the large- Ai'c suppressions to a certain extent. We perform the fit of the L[ using the new 
set C^^^^"^^ and we check the obtained. If the decreases then we substitute the Cl^°^^^ 
with Sometimes we also allow (j^^^™^ to be selected even though the corresponding 

is not smaller than the previous one (see last step in the flow diagram of Figured]). This 
is done so to let the to take quite different values and thus to test as many different sets 
as possible. It is also needed for our algorithm to be able to move out of local minima. We 
find that, when we let our algorithm run long enough, we cover quite many different sets 
of C[. We chose different initial to widen their range of variability. In addition we have 
started the random walk from the same starting point several times including different 
random starting points. We chose as random starting points l/3/(167r^)^ with the extra 
factor of 1/3 since without the extra 1/3 we never reached a smaller than one. 

The fits are performed including the same input as for fit All but with a few extra 
requirements. We add as input the curvature of the scalar formfactor in fl30]) and all 
the ii of the last column of Table El We do not instead demand to fit all the tttt and ttK 
scattering parameters, since it costs in terms of computing time. So, as done for fit All, we 
include only Oq, a^, al^'^ and a^^'^. Notice that we do not find any large discrepancies when 
we add more scattering parameters in fit All, as was remarked at the end of Section El 

We also require a good convergence of the masses and decay constants expansions. 
The reason is that in this way we have the possibility to "select" those granting us 
convergence for those quantities. This also allows to keep under control the quality of the 
fits. Otherwise too much freedom would be left to the constants, and many different fits 
with a low but with very bad convergences, can be reached. Clearly such convergence 
constraints have a strong effect on the L[ constants as described in the next Section 17.11 

7.1 Convergence constraints 

We devote this section to a discussion of the convergence constraints imposed on the masses 
and decay constants in the fits with random C^. Let us first show the case of the decay 
constants. 
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We performed the fits constraining the NNLO contributions to Ft^, Fk and Fk/Ft^ 
constants to be small, i.e. less than the 10% of the LO ones. Remember that the expansions 
for the and Fk decay constants are 

Fk 
F^ 



Fo+FkL^ + F, 



K\ri6 , 



and that in our fits, as explained in Section 13.21 we include their ratio as 



(43) 



1 + 



A -TO 



Fo 



Fk 
Fo 



F^ 
Fo 



,4 -To 



NLO 



NNLO 



(44) 



Specifically the convergence constraints are included through the following partial 



i(part) ■ 



Xi{part) 



Xj(part) 




/0.05 



^ /0.07 

NNLO 



Xi(part) 




/0.05 



(45) 



With ( l45l) we have been able to find many fits with a good convergence. On the other hand 
the NLO contribution for F^^ turns out to be smaller than the expected 30%. The reason 
resides in the third of the relations in ( H5l) . Requiring all the NNLO pieces for all the decay 
constants to be small implies that the single contributions Fj^lpe and Fk\p6 are small. But 
also the term (F^/Fo)|p4 must be small, otherwise the NNLO contribution of Fk/Fj^ is 
allowed to be large. This leads to small NLO corrections for F^^ and thus a Fq ~ F^^. 
We apply similar restrictions also to the masses 



Xi(part) 



m 



M 



m 



MO 



/O.l 



(46) 



where M stands for vr, K and t] mesons and mMo are the leading order contributions to 
the masses. We have kept the value of ms/m fixed at 27.8 as for fit All. 

We conclude this section with a final remark. One might wonder why we have not 
imposed similar constraints also for fit All, since these could improve the convergence of 
the expansions. The reason is that when we require them, we obtain a reasonably good fit 
[x^ 8 with 10 degrees of freedom ) and with better constrained L4 and Lg. But it also 
causes much worse predictions for the e.g. £3 ~ 6.5. 



7.2 Results 

Now we are ready to show the outcomes of our studies when the Cl are set to random 
values using the procedure of Figure [Hand with the constraints listed in Section [711 above. 
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^r(in) 


resonance 


zero 




random 


1 7'/ 


yv.v, i.z,) 


('0 5 1 1 


(0 9 Of,] 


fO 3 19) 




(0.4,1.2) 


(0.2,0.8) 


(-0.2,0.5) 


(0,1.6) 




(-5.5,-3.0) 


(-4.5,-3.0) 


(-3.3,-1.2) 


(-6.2,-2.0) 




(0.15,0.35) 


(-0.1,0.25) 


(-0.1,0.2) 


(0.05,0.35) 




(1.2,1.6) 


(1.35,1.5) 


(1.25,1.55) 


(1.2,1.6) 




(-0.05,0.25) 


(-0.2,0.2) 


(-0.15,0.15) 


(-0.05,0.3) 




(-0.45,-0.1) 


(-0.43,-0.18) 


(-0.43,-0.24) 


(-0.45-0.2) 




(0.4,0.8) 


(0.45,0.72) 


(0.5,0.731) 


(0.45,75) 



Table 15: The ranges for L[ values as obtained changing the CI according to a random 
walk algorithm. The different ranges correspond to different initial values for the C^. With 
CQM we indicate the CI of the chiral quark model. All the fits have a < 1- For all the 
fits /i = 0.77 GeV. See the description in the text for further details on how the fits have 
been performed. 



First of all we must point out that due to the freedom we allow to the many different 
fits of the L[ have a low x^. We set initial CI equal to zero, resonance exchange or chiral 
quark model estimates as well more random starts. We easily reach < and we found 
many fits with x^ ^ 0.5. Reducing the steps of the random walk we can even find smaller 
values. However once the x^ reaches a reasonably low value, e.g. 1, there is no apparent 
reason why one should prefer one specific fit to another. Due to the several different sets 
of CI under study, we can only quote the ranges where the L[ vary and where we obtain 
a < 1- Such ranges are quoted in Table [15] and are obtained, starting from different 
initial sets Cj''™\ Keep in mind that these ranges depend on the C[ chosen. Plugging 
in a LI fit without the corresponding C[ will not produce any sensible results. The way 
we determined those numbers is shown in Figure [2] on the example for L\ where we have 
plotted a number of fits that gave x^ <1 for the different starting points. We have typically 
stopped the fits when a x^ of about 0.4 or below was found and the tails at low x^ are an 
artefact, they were done with runs with a very low e and a very low T. 

The results are rather cumbersome. Lg look quite free to vary in large 

ranges. Notice also that the large relation 2L\ ^ L2 is still not recovered. On the 
other hand, due to the converge constraints of Section 17.11 we narrow the intervals for the 
other constants. Especially L5 takes a large value and L\ a small one. As explained in 
Section \7~A\ we essentially require that Fk/Ft^, Fx and F^^ have small NNLO corrections 
which give that F^/ F^^ is given by the NLO and thus determines L5 at a fairly large value. 
That has small corrections at all then in turn requires a fairly small L\. The choice 
to constrain the convergence of those quantities is dictated by the lack of information to 
constrain more the CI- When we release such constraints indeed we find different looking 
fits, but affected by a bad convergence. Somewhat more surprising is that the Lg typically 
takes on values that are smaller than Lg. 

As you can see no clear final conclusion can be drawn with such results. When we 
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Figure 2: The values of L\ for the random walk fits of Table [151 In the plot all the fits are 
collected. The ranges of variability for is quite large. In the picture it is also evident 
how the two values of the couplings depend somwehat on the different initial values of the 
constants. The are at the scale /i = 0.77 GeV. 
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best reso 


best random 






n sc: -1- n DQ 
u.otj HI u.uy 










-3.91 ±0.28 


-3.51 ±0.28 




0.16 ±0.10 


0.20 ±0.10 




1.40 ±0.09 


1.40 ±0.09 




0.10 ±0.14 


0.12 ±0.14 




-0.32 ±0.13 


-0.32 ±0.13 




0.64 ±0.16 


0.63 ±0.16 




0.30 


0.36 



Table 16: The results as obtained using the CI from the best found starting from the 
resonance or from a completely random one as described in the text. The L\ are given at 
/i = 0.77 GeV. The corresponding CI sets can be found in Table [18] in appendix. 



d 


best reso 


best random 






4 

p 






4 

P 




ml 


0.987 


0.021 


-0.008 


0.993 


0.021 


-0.012 


m\ 


1.057 


-0.054 


-0.003 


1.060 


-0.058 


-0.002 


m\ 


1.132 


-0.133 


0.001 


1.136 


-0.135 


-0.001 


1 


0.178 


-0.010 


1 


0.187 


-0.010 


FrI Fq 


1 


0.395 


0.009 


1 


0.404 


0.011 


Fk/F^ 


1 


0.217 


-0.020 


1 


0.217 


-0.020 



Table 17: The convergence for the best found starting from the resonance estimate, 
here Fq = 0.079 GeV. Fit as in left column of Table [161 The convergence for the best 
found starting from the fully random estimate, here Fq = 0.078 GeV. Fit as in left column 
of Table M 

performed this study we were hoping not to find as many good fits and smaller ranges for 
the L[. The study shows instead that it is very difficult, if not impossible, to narrow the 
ranges for the L[ with such a poor knowledge of the C[. On the other hand it also shows 
that fits of the L[ with good convergence do exist, if the C[ are changed. In Table [T6l we 
show the L[ obtained for the smallest found starting from the resonance estimate and 
from fully random C[ as described above and we show the convergence for some quantities 
in those two fits in Table [T71 The fits with very low we have obtained, such as the two 
shown here, tend to have similar expansions for the masses and the decay constants. In 
order to see how the various CI' look like we have added the values for these two fits in the 
appendix. 

We can draw some conclusions by studying correlations. The effect of £2 is very visible if 
we plot for the various fits with < 1 L2, which is the NLO dependence of £2 on the LECs, 
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1 .6e-05 




Figure 3: The correlations between and 2C[3 — Cn- This resuhs from including the 
value of £2 in the fit. 

versus 2Cl^ — C[i, which is the dependence on the NNLO LECs. The tight correlation 
shows that the NNLO contribution depends very little on the other L[. More precisely, 
this shows the NLO LEG combination at NLO versus the NNLO LEG combination at 
NNLO and that the NNLO contribution depends fairly little on the value of the L[. We 
get more of those constraints directly. A very similar one is from the value of ii shown 
in Figure |H The correlations for other observables tend to be weaker indicating that the 
NNLO contributions are more dependent on the value of the L[ for these cases. We show 
examples with a weaker but still existing correlation in Figure |5] where the correlations 
resulting from F^/Fj^ are shown and in Figure IH] for (r^)^. In both cases we have plotted 
on the horizontal axis the combination of L[ the quantity depends on at NLO and on the 
vertical axis the combination of C[ the quantity depends on at NNLO. 

There are also correlation between the fitted values of the L[. L\, and L3 show a 
reasonable correlation among themselves but are essentially not correlated with the others. 
There are weaker correlations between L4 and Lg and between Ly and Lg. These correla- 
tions are shown in Figures [71 El El CHI and [TTJ We have shown a curve in all plots guiding 
the eye as well and given it in the caption of the figure. 

Note that throughout this section we have considered all fits with a < 1 to be 
essentially possible. 
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2.5e-05 




-1e-05 I 1 1 1 1 1 1 1 

-0.004 -0.0035 -0.003 -0.0025 -0.002 -0.0015 -0.001 -0.0005 



Figure 4: The correlations between 2L\ + Lg and Cg — Cli + 4C[3. This results from 
including the value of ii in the fit. 



-7e-06 




-1.5e-05 I 1 1 1 1 1 1 1 1 1 1 

0.00115 0.0012 0.00125 0.0013 0.00135 0.0014 0.00145 0.0015 0.00155 0.0016 0.00165 



Figure 5: The correlations resulting from Fk/F^ 
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Figure 7: The correlations between L[ and L2, the curve shown is L2 = 1-4(L^ — 0.00032). 
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-0.001 




-0.0065 I 1 1 1 1 1 1 1 1 1 1 

0.0002 0.0003 0.0004 0.0005 0.0006 0.0007 0.0008 0.0009 0.001 0.0011 0.0012 

Figure 8: The correlations between L'[ and Lg, the curve shown is Lg = —4.7(1/^ + 0.00003). 



-0.001 
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-0.0002 0.0002 0.0004 0.0006 0.0008 0.001 0.0012 0.0014 0.0016 0.0018 



Figure 9: The correlations between and Lg, the curve shown is Lg = — 3. 1(L2 + 0.00055). 
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Figure 10: The correlations between L4 and Lg, the curve shown is Lg = L4 — 0.0001. 

0.0008 I 1 1 1 1 1 1 1 1 1 




0.00035 I 1 1 1 1 1 1 1 1 1 

-0.0005 -0.00045 -0.0004 -0.00035 -0.0003 -0.00025 -0.0002 -0.00015 -0.0001 -5e-05 



Figure 11: The correlations between and Lg, the curve shown is Lg = —{Lj — 0.0003). 
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8 Conclusions 



In this paper we have shown the results for a new global fit of the L[ at NNLO, with 
techniques similar to the ones in [HI |9] . Different treatments of the coupling constants 
have been considered: the resonance estimate of [8], the results of [13] and the use of 
randomly selected C[. The results are difficult to interpret and unexpected. All the fits 
that have been performed using the NNLO couplings from [8] or [13] show both strong and 
weak points. The fits obtained from the randomly selected C[ are too different from each 
other to draw a final conclusion, although they give a rough indication on where we can 
expect to find the values of L[. They also provide a proof of principle that with reasonable 
values of the C[ a reasonably convergent series for SU (3) ChPT can be obtained for many 
quantities. 

The fit that presents the least discrepancies and best convergence of the chiral expansion 
is fit All in Table El which has been obtained with the resonance estimate of the C^. It 
succeeds in fitting many observables like the vrvr and nK scattering parameters and the 
slope of the scalar formfactor of the pion. It also reproduces quite well the experimental 
results for the fs and gp K^^^ formfactors although it does not predict the curvature of fs- 
The perturbative expansions for masses and decay constants reported in (l36ll and (1371) look 
suspicious but are acceptable. On the other hand it does not satisfy the large N^. relation 
2L\ and it fails in well constraining the and Lg values. Finally its prediction for 

^2 is far from the current estimate of that constant. We have at present not included many 
results from lattice QCD. We expect that this should improve in a few years allowing for 
another step forward in confronting ChPT with data. 
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A Q values 

The aim of this section is to present the values used for the C[. We only present results 
for the CI that actually contribute to the observables we have included in the fits and we 
indicate with the superscript * the ones that are subleading in N^.. They can all be found 
in Table UHl The column labelled reso is the resonance exchange estimate of Section HI 
The column labelled CQM is the estimate of [13j as discussed in Section 16. 4[ The values 
directly from their model can be found in Table IV of The values in the column 

labelled CQM have been multiplied with nomalization factor a = 0.27 from the fit in 
the right column of Table l9l The normalization from the fit to the linear NA48/2 input. 
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Table [TT] is a little smaller but essentially the same. The last two columns are from the 
random walk/simulated annealing estimates for the CI reported in Section [7] from the best 

found starting from the resonance estimate and a fully random starting point. These 
are the CI for the fits reported in Table [161 

The main purpose is to show the typical sizes of the C[ we obtained for the fits and 
that the pattern for the good fits can be quite different. 

References 

S. Weinberg, Physica A96 (1979) 327. 

J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 

J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 

J. Bijnens, G. Colangelo, G. Ecker, JHEP 9902 (1999) 020. [hep-ph/9902437]. 

J. Bijnens, Prog. Part. Nucl. Phys. 58 (2007) 521 [hep-ph/0604043]. 

J. Bijnens, I. Jemos, Eur. Phys. J. C64 (2009) 273-282. [arXiv:0906.3118 [hep-ph]]. 

J. Bijnens, G. Colangelo, J. Gasser, Nucl. Phys. B427 (1994) 427-454. [hep- 
ph/9403390]. 

G. Amoros, J. Bijnens and P. Talavera, Nucl. Phys. B 585 (2000) 293 [Erratum-ibid. 
B 598 (2001) 665] [hep-ph/0003258]. 

G. Amoros, J. Bijnens and P. Talavera, Nucl. Phys. B 602 (2001) 87 [hep-ph/0101127]. 

J. Bijnens, P. Dhonte, JHEP 0310 (2003) 061. [hep-ph/0307044]. 

J. Bijnens, P. Dhonte, P. Talavera, JHEP 0401 (2004) 050. [hep-ph/0401039]. 

J. Bijnens, P. Dhonte and P. Talavera, J. High Energy Phys. 0405 (2004) 036 [hep- 
ph/0404150]. 

S. -Z. Jiang, Y. Zhang, C. Li et al, Phys. Rev. D81 (2010) 014001. [arXiv:0907.5229 
[hep-ph]]. 

G. Ecker, P. Masjuan, H. Neufeld, Phys. Lett. B692 (2010) 184-188. [arXiv: 1004.3422 
[hep-ph]]. 

A. Pich, Lectures at Les Houches Summer School in Theoretical Physics, Session 68: 
Probing the Standard Model of Particle Interactions, Les Houches, France, 28 Jul - 5 
Sep 1997, [hep-ph/9806303]. 

S. Scherer, Adv. Nucl. Phys., 27 (2002) 277 [hep-ph/0210398]. 



37 



i 


reso 


CQM 


best 


best 


i 


reso 


CQM 


best 


best 








reso 


rand 








reso 


rand 


1 


1.216 


0.866 


1.683 


-0.733 


25 


-1.838 


-1.366 


-1.452 


1.282 


2* 


0.000 


0.000 


0.113 


0.280 


26 


-0.284 


0.765 


-0.597 


-0.485 


3 


0.000 


-0.011 


0.324 


0.084 


27* 


-0.261 


-0.352 


-0.228 


-0.155 


4 


1.452 


0.708 


2.225 


1.266 


28* 


0.135 


0.069 


0.177 


0.147 


5 


0.619 


-0.231 


0.779 


1.147 


29 


-1.363 


-0.704 


-1.907 


-0.785 


6* 


0.000 


0.000 


-0.307 


-0.050 


30* 


0.270 


0.137 


0.165 


0.545 




0.000 


0.000 


0.350 


-0.003 


31 


-0.616 


-0.144 


-0.389 


1.310 


8 


0.619 


0.528 


1.434 


0.615 


32* 


-0.002 


0.041 


0.291 


0.356 


9* 


0.000 


0.000 


0.148 


-0.194 


33* 


0.208 


0.021 


0.291 


-0.102 


10 


-1.239 


-0.240 


-0.164 


0.307 


34 


1.432 


0.363 


2.321 


1.077 


11* 


0.000 


0.000 


-0.112 


-0.344 


35* 


-0.009 


0.039 


0.236 


0.146 


12 


-0.619 


-0.078 


-1.358 


-0.512 


36* 


0.000 


0.000 


-0.169 


-0.087 


13* 


0.000 


0.000 


0.265 


-0.002 


63 


0.619 


0.683 


0.665 


0.776 


14 


0.000 


-0.190 


-0.759 


-0.828 


64* 


0.000 


0.000 


0.389 


-0.095 


15* 


0.000 


0.000 


-0.228 


-0.233 


65 


1.239 


-0.555 


0.611 


0.432 


16* 


0.000 


0.000 


0.007 


0.063 


66 


1.049 


0.391 


1.703 


0.416 


17 


0.000 


0.002 


0.125 


1.121 


67* 


0.000 


0.000 


-0.304 


0.017 


18* 


-0.202 


-0.128 


-0.284 


-0.063 


68* 


0.000 


0.000 


0.002 


0.536 


19 


0.001 


-0.110 


-0.411 


-1.147 


69 


-0.577 


-0.196 


-0.664 


-0.784 


20* 


-0.002 


0.041 


-0.335 


-0.043 


83 


0.163 


0.016 


-0.294 


-0.553 


21* 


0.001 


-0.014 


0.018 


-0.088 


84* 


0.000 


0.000 


0.357 


-0.264 


22 


-0.297 


0.062 


0.545 


1.117 


88 


-1.383 


-1.249 


-0.912 


-0.331 


23* 


0.000 


0.000 


0.048 


0.269 


90 


5.069 


0.557 


5.238 


-0.204 


24* 


0.811 


0.370 


0.694 


0.013 













Table 18: The C[ constants as obtained by several methods. In the table we quote 10^ x C[ 
were the i is given in the first column. The i* indicates a Nc suppressed C^. The second 
column corresponds to the resonance model of Section HI In the third one we quote the 
of [13] multiplied by a = 0.27. The particular value of a is chosen fitting it to the input 
observables together with the L[, as explained in Section 1^31 In the last two columns we 
quote two sets of C[ corresponding to the best fits of the L[ when the are realesed, as 
explained in Section [71 The column labeled best reso comes from initial values of the 
equal to the ones of the resonance estimate, while in the column labeled best rand even 
the initial values of the are chosen randomly. The scale of renormalization for all the 
sets is /i = 0.77 GeV. 
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